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Abstract: We evaluate neutralino-nucleon scattering rates in several well- motivated su- 
persymmetric models, and compare against constraints on the neutralino relic density, 
BF(b — > 57) as well as the muon anomalous magnetic moment a^. In the mSUGRA 
model, the indirect constraints favor the hyperbolic branch/focus point (HB/FP) region of 
parameter space, and in fact this region is just where neutralino-nucleon scattering rates 
are high enough to be detected in direct dark matter search experiments! In Yukawa unified 
SUSY 50(10) models with scalar mass non- universality, the relic density of neutralinos is 
almost always above experimental bounds, while the corresponding direct detection rates 
are below experimental levels. Conversely, in five dimensional SO (10) models where gauge 
symmetry breaking is the result of compactification of the extra dimension, and super- 
symmetry breaking is communicated via gaugino mediation, the relic density is quite low, 
while direct detection rates can be substantial. 
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1. Introduction 



A wide variety of astro-physical data now point conclusively to the existence of cold dark 
matter (CDM) in the universe. The most recent results come from the Wilkinson Microwave 
Anisotropy Probe (WMAP)Q. Their results confirm the standard model of cosmology and 
fit its parameters to high precision. The properties of a flat universe in the KCDM model 
are characterized by the density of baryons matter density (fi m ), vacuum energy (Oa.) 
and the expansion rate (h) which are measured to be: 

Q b = 0.044 ± 0.004 (1.1) 
tt m = 0.27 ±0.04 (1.2) 
n A = 0.73 ±0.04 (1.3) 

h = 0.71+8:8* (i.4) 

From the WMAP results one derives the following value for the cold dark matter density: 

n CDM h 2 = o.ii26i8;889o(±8:8i8i) a * 68(95)% cl. (1.5) 

In spite of the excellent fit to the data, the standard cosmological model has several 
fundamental open questions. One of these is the origin and nature of dark matter. Within 
the context of imparity conserving supersymmetry (SUSY) j| the lightest supersymmetric 
particle (LSP) offers a robust solution to this problem. First, SUSY itself is the most 
complete theoretical extension of the Standard Model (SM), and is well motivated experi- 
mentally. This solidly supports the neutralino LSP dark matter theoretically. Second, the 
WMAP results strongly restrict neutrino hot dark matter to Q u h 2 < 0.0076 (95%CL.) and 
also exclude warm dark matter based on detection of re- ionization at z ~ 20. Warm dark 
matter in the form of gravitinos could also be a good candidate for dark matter in certain 
supersymmetric models, but WMAP has vetoed this option. Therefore, in the light of 
recent WMAP results, the neutralino LSP emerges as one of the most attractive candidate 
for CDM. 

If indeed all space is filled with relic neutralinos, then it may be possible to directly de- 
tect them via their scattering from nuclei, or to indirectly detect them via their annihilation 
products. There are two kinds of non-accelerator experiments aimed at the CDM search. 
One way is to search for the products of neutralino annihilation in outer space or inside the 
sun or earth where they become concentrated due to the gravitational force. Another way 
to search for CDM is via direct detection from their scattering off nuclei by measuring the 
nuclear recoil. There are several existing and future projects engaged in the direct search 
for weakly interacting massive particles (WIMPS) [||. These include experiments based 
on germanium detectors such as IGEX §, HDMS §, CDMS §, EDELWEISS and 
GENIUS §]. Scintillator detectors are used by DAMA [§ and ZEPLIN 0, 0, @ 0. 



Projection chambers utilized in DRIFT jl4|], and metastable particle detectors in SIM- 
PLE @ and PICASSO @. 



Presently the best limits on the spin-independent WIMP-nucleon cross-section (crsi) 
have been obtained by the CDMS, EDELWEISS and ZEPLIN1 groups, while a signal is 



- 1 - 



claimed by the DAMA collaboration. Collectively, we will refer to the reach from these 
groups as the "Stage 1" dark matter search. Depending on the neutralino mass, the com- 
bined limit on the neutralino-proton spin-independent cross section asi varies from 10 -5 
to 1CT 6 pb. This cross section range is beyond the predicted levels from most supersym- 
metric models if one also requires the model to be within the experimental constraints from 
LEP2, (g - 2)^, BF(b — ► 57), relic density QcDAih 2 and BF(B S — > [i + data. How- 
ever, experiments in the near future like CDMS2, CRESST2, ZEPLIN2 and EDELWEISS2 
(Stage 2 detectors) should have a reach of the order of 10 -8 pb. Finally, a number of 
experiments such as GENIUS, ZEPLIN4 and CRYO ARRAY are in the planning stage. We 
refer to them as Stage 3 detectors, which promise strong limits of the order of asi < 10 -9 
- 1CP 10 pb, and would allow the exploration of a considerable part of parameter space of 
many supersymmetric models. 



Beginning with the paper by Goodman and Witten |17], there have been numerous 



studies on neutralino nuclei scattering rate evaluation [18, |l^, 2C, 21, 22, E^], with the 
trend of gradually improving the quality of calculations and extending the range of super- 
symmetric models. It is useful to note that in the case of the general MSSM, asi could 
be several orders of magnitude higher |2^] than in more constrained models such as the 
minimal supergravity model (mSUGRA). 

In this paper, we perform an updated analysis of neutralino elastic scattering off nuclei 
in a class of supersymmetric models. There are various motivations for this study. 

• We have included constraints from the neutralino relic density the rare decays 

b — > sj and B s — > fi + fi~, the muon anomalous magnetic moment a^, as well as from 
the LEP2 experiment. Together these constraints significantly restrict the SUSY 
model parameter space in the relevant range of asi- 



We are using a new version of ISAJET, version 7.65[24], in which the complete one- 
loop corrections to sfermion masses |2^] have been incorporated. In addition, the 
stability of the new version has been greatly improved in the low region which 
allows us to better access the hyperbolic branch/focus point (HB/FP) region of the 
mSUGRA model. This is important, because the HB/FP region can simultaneously 
satisfy constraints from fij^/i 2 , (g — 2) M , and BF(b — > sj) while offering a decoupling 



solution to the SUSY flavor and CP problems, and maintaining naturalness [26, 27 



Besides the mSUGRA model, we also analyse SUSY grand unified theories (GUTs) 



based on the 50(10) group [29] where the MSSM Higgs doublets are both present 
in the same 10-dimensional Higgs multiplet. These models predict unification of the 
Yukawa couplings for the third generation: ft = ft = fr- Following the strategy of 



our previous study [3C], we examined models with a high degree of Yukawa coupling 
unification, and have evaluated neutralino nuclei scattering rates together with ex- 
perimental constraints. For Yukawa unified models with fi < 0, there exist parameter 
choices consistent with neutralino relic density constraints, although direct detection 
rates are frequently quite low. For models with [i > 0, the relic density is almost 
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always above measured limits, while direct detection rates are almost always below 
experimentally accessible levels. 



We also examine an SO(W) motivated SUSY GUT model [31] with non-universal 



gaugino masses and otherwise vanishing soft-parameters at the GUT scale. This 



setup is well motivated by SUSY GUTs constructed in higher dimensions [32, 33] 



utilizing gaugino mediated SUSY breaking. Many of these models elegantly solve the 



doublet-triplet splitting and fast proton decay problems of 4D SUSYGUT models [34] 
by the oribifold compactification of the extra dimenstions. In the phenomenologically 
viable region of the parameter space of this model the lightest neutralino generally 
has a large higgsino or wino component, which leads to low values of il^^h 2 , and 
neutralino dark matter would need to be augmented by other forms of CDM, such as 
axions, or hidden sector states. We find that a si can frequently be several orders of 
magnitude larger than the mSUGRA case, owing to the large wino and/or higgsino 
component of the neutralino. 

The rest of this paper is organized as follows. In Section 2 we discuss details of 
our evaluation of the neutralino-nucleon elastic scattering cross section. In Section 3, we 
present results for the mSUGRA model, and point out the intriguing features endemic 
to the HB/FP region. In Section 4, we perform studies for 50(10) motivated Yukawa 
unified SUSY GUT models, while Section 5 presents results for SUSY GUT models with 
non- universal gaugino masses. In Section 6 we present our conclusions. 

2. Details of calculation for neutralino elastic scattering on nuclei 

The interactions for elastic scattering of neutralinos on nuclei can be described by the sum 
of spin- independent {C e J c l lm ) and spin-dependent {C e J p ( n ) Lagrangian terms: 

£ e ff. =C eff , +C eff . (2.1) 

elastic scalar ' spin V / 

In this paper we evaluate the spin-independent cross section a$i of neutralino scattering off 
of nuclei which is the main experimental observable since asi contributions from individual 
nucleons in the nucleus add coherently and can be expressed via SI nuclear form-factors. 
The cross section asi receives contributions from neutralino-quark interactions via squark, 
Z and Higgs boson exchanges, and from neutralino-gluon interactions involving quarks, 
squarks and Higgs bosons at the 1-loop level. The differential asi off a nucleus X% with 
mass rriA takes the form pq] 



da SI _ _L_ 
d\q\ 2 irv 2 



[Zf p +(A-Z)f n ] 2 F 2 (Q r ), (2.2) 



where q = — +m l v is the three-momentum transfer, Q r = and F (Q r ) is the scalar 

z i 

nuclear form factor, v is the velocity of the incident neutralino and f p and f n are effec- 
tive neutralino couplings to protons and neutrons respectively. This formalism has been 
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reviewed in [35, [H], |2(|. Explicit expressions for f p and f n can be found, e.g. in |2C]. The 



original calculation has been done in [19] and can be expressed as 



m N 



AN) 

JTq 



m„ 



q=u,d,s 



-m q m^g q 



_f_ AN) \- 

+ 97 JTG 



27 



c,b,t 



+ 



(2.3) 



where N = p, n for neutron, proton respectively, and = 1 — Ylq=u d s Itq ■ The 

expressions for the f q couplings as well as other terms denoted by • • • are omitted for the 



sake of brevity but can be found in [19 



(p) 

The parameters frf q , defined by 



< N\m q qq\N >= m N f^ (q = u,d,s) 



(2.4) 



contains uncertainties due to errors on the experimental measurements of quark masses. 

(p) 

We have adopted values of renormalization-invariant constants fj,' and their uncertainties 
determined in [21] 



= 0.020 ± 0.004, J® = 0.026 ± 0.005, fj$ = 0.118 ± 0.062 (2.5) 



(n) 



Tu 



0.014 ± 0.003, 



(n) 



Td 



0.036 ± 0.008, 



(n) 
Ts 



0.118 ±0.062. 



(2.6) 



In this paper we calculate the quantity which is being conventionally used to compare 
experimental and theoretical results - the cross section <jp for neutralino scattering off 
the proton in the limit of zero momentum transfer 



-m N2 f 2 
rn r Jn 

TT 



(2.7) 



where m N = m^m^J(vn jy + Tn^). 

In our calculations we have used the CTEQ5L set of parton density functions p6| 
evaluated at the QCD scale Q = yM^ USY — m 2 ~ . 



3. Direct detection rates in the mSUGRA model 



We start with the analysis of the neutralino scattering rates in the minimal supergravity 
model [37]. The mSUGRA model assumes universal boundary conditions at the GUT scale 
for scalar and gaugino masses as well as for trilinear ^4-parameters and therefore is defined 
by a set of just four parameters and a sign: 



m 0i m i/2i A)> tan /3 and sign (/i), 



(3.1) 



where tan/3 = v u /vd parametrizes the Higgs sector. The top-quark pole mass is taken to 
be m t = 175 GeV. 
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3.1 Experimental constraints 

Before presenting results we discuss the following experimental constraints on the mSUGRA 
model. 



LEP2 constraints. Based on negative searches for superpartners at LEP2 38, 
we require > 103.5 GeV and R > 99 GeV provided m^ — m^ i > 10 GeV] 
which is the most stringent of the slepton mass limits. The LEP2 experiments also 
set a limit on the SM Higgs boson mass: fnu SM > 114.1 GeV f40f| . This is of relevance 
since in our mSUGRA parameter space scans, the lightest SUSY Higgs boson h is 
almost always SM- like. We implement the MSSM Higgs boson mass bounds from 
Ref. [fil[], although these bounds co-incide with the SM bound since in our mSUGRA 
scans niA > 120 GeV. 

Neutralino relic density. The WMAP results Q set a stringent bound on the 
neutralino relic density, given by 0.0945 < QcDAih 2 < 0.129 at 95% CL. The upper 
limit above represents a true constraint, while the corresponding lower limit is flexible, 
since there may be additional sources of CDM such as axions, or states associated 
with the hidden sector and/or extra dimensions. To estimate the relic density of 
neutralinos in the mSUGRA model, we use the recent calculation in Ref. |fffi| . In that 
work, all relevant neutralino annihilation and co-annihilation reactions are evaluated 
at tree level using the CompHEP |43[| program. The annihilation cross section times 
velocity is relativistically thermally aver aged [|44|1, which is important for obtaining 
the correct neutralino relic density in the vicinity of annihilations through s-channel 
resonances. 

b — ► S7 decay. The branching fraction BF{b — > s~f) has been measured by the 
BELLE Q, CLEO[g6) and ALEPHQ collaborations. A weighted averaging of re- 
sults yields BF(b -» s 7 ) = (3.25 ± 0.37) x 10" 4 at 95% CL. To this we should add 
uncertainty in the theoretical evaluation, which within the SM dominantly comes 
from the scale uncertainty, and is about 10%. Together, these imply the bounds, 
2.16 x 10~ 4 < BF(b — ► sj) < 4.34 x 10 -4 . In our study, we show contours of 
BF{b —* sj) allowing the reader to decide the extent to which parameter space is 
excluded fig] . The calculation of BF(b — * 57) used here is based upon the program 
of Ref. [j49f| . Our value of the SM b — > sj branching fraction yields 3.4 x 10~ 4 , with 
a scale uncertainty of 10%. 

Muon anomalous magnetic moment. The muon anomalous magnetic moment 
a n = (9 ~ 2)^/2 has been measured to high precision by the E821 experiment p0|. 
Comparison of the measured value against theoretical predictions gives, according 
to Hagiwara et al. [|5l]] : 11.5 < 5a^ x 10 10 < 60.7. A different assessment of the 



theoretical uncertainties [51] using the procedure described in Ref.[|4^] gives, —16.7 < 



5a a x 10 < 49.1. In view of the theoretical uncertainty, we only present contours 
of 8a n, as calculated using the program developed in [p^j , and leave it to the reader 
to decide the extent of the parameter region allowed by the data. 
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B s —* u + u decay. The branching fraction of B s to a pair of muons has been 
experimentally bounded by CDF|j|: BF(B S -> n + n~) < 2.6 x 1CT 6 . While this 



branching fraction is very small within the SM (BFsm(B s —* u + u ) ~ 3.4 x 10 9 ), 
the amplitude for the Higgs-mediated decay of B s grows as tan 3 j3 within the SUSY 
framework^], and hence can completely dominate the SM contribution if tan/3 is 
large. In our analysis we use the results from [55] to delineate the region of mSUGRA 



parameters excluded by the CDF upper limit on its branching fraction. 
3.2 mSUGRA results 

In Fig. El we show our first results of the spin independent neutralino-nucleon elastic 
scattering cross section in the m,Q vs. m^ 2 parameter plane for tan/3 = 10 — 55, with 
Aq = and for u < (upper two rows) and /i > (lower two rows). 

In this figure we have applied only LEP2 constraints (black color) and theoretical 
constraints denoted by red color which show the regions forbidden either due to lack of 
REWSB (lower right region) or because f\ is the LSP (upper left region). 

One can see that for both positive and negative /i the the highest value for asi is about 
10~ 6 pb. This happens in two regions: either in the lower left corner of the tuq vs. m\/2 
parameter space where the exchanged squark masses are light, or in the HB/FP region for 
large mo values and low-intermediate values (gray and green colors), along the border 
of the red forbidden region where REWSB does not occur. In this region, there is a large 
higgsino component to the neutralino, which enhances the Higgs exchange contribution 
to o"s/^6|. For u > 0, the lowest cross sections reach below 10~ 10 pb and occur at very 
large values of uiq and m 1 / 2 where squark masses are high, thus suppressing the overall 
scattering cross section. For \x < 0, the lowest cross section values can dip far below 
10~ 12 pb in the magenta-shaded regions. The nature of this behavior in asi is due to 
a cancellation between processes involving up- and down-type quarks where the leading 
contribution comes from i-channel Higgs boson exchange pl|]. 

In Fig. |2| we plot asi versus neutralino mass for a random scan over mSUGRA pa- 
rameter space for various tan [3 and Aq = 0, for /x < (top frame) and \i > bottom 
frame). Parameter space points which satisfies only theoretical and LEP2 constraints (i.e. 
parameter space presented in Fig. p]) are represented by small dots. If instead we restrict 
&z h 2 to be within the experimental limits 0.0945 < QcDAih 2 < 0.129, only the solid 
circles survive. The different colors of the dots correspond to different values of tan (3. We 
can see now how restrictive the relic density constraint is for the parameter space and for 
the range of asi- One can see that the maximum value asi can be reduced almost by an 
order of magnitude for some values of tan (3. We also show in these plots the reach of Stage 
1, Stage 2 and Stage 3 experiments. Evidently Stage 1 experiments are just now starting 
to explore regions of mSUGRA model parameter space 1 57]. Stage 2 and Stage 3 detectors 



can explore much larger regions of parameter space, but none of the planned detectors will 
be able to completely rule out SUSY dark matter. 

The regions of parameter space allowed by the relic density constraint form branches 
of distinct patterns. One sort of branch is the HB/FP region. It is exhibited by region of 
almost constant asi ~ 10~ 8 pb which is nearly independent of tan (3. In this region, it is 
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Figure 1: Cross section rates for spin-independent neutralino-proton scattering in the mSUGRA 
model for Aq = 0. 



possible to satisfy not only the relic density constraints, but also (as shown below) and 
BF(b — ► S7) constraints. One can solve the SUSY flavor and CP problems, and possibly 
maintain naturalness. The Stage 3 detectors ought to have sufficient reach to cover this 
interesting region of model parameter space! It is important to note, however, that in 
the HB/FP the relic density can be quite low, so that the assumed local density of dark 
matter should be rescaled[p^]. Hence, in all subsequent figures we will plot / • a si, where 
/ = O^/i 2 /0.094 for n^h 2 < 0.094, and / = 1 for Vt^h 2 > 0.094. Our reach plots will 
also be determined in terms of / • a si rather than in terms of a si- 

In Fig. ||a.), we show / • a si vs. as a function of the complete model parameter 
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Figure 2: (75/ versus neutralino mass 



space for \i < 0. The ranges of parameters are listed above the figure, and include a scan 
over all Aq values. We also show the reach of Stage 1, 2 and 3 detectors. Most points 
in mSUGRA model parameter space are now excluded by the WMAP constraint, so we 
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m z{ (GeV) Or 

• Oh 2 < 0.094 • 0.094 < Qh 2 < 0. 129 • 0. 129 < QJi 2 
LEP2 Stage 1 ■ - Stage 2 — Stage 3 

Figure 3: In a) we plot rescaled cross section rates vs. for a scan over mSUGRA parameter 
space with fi < 0. In 6), we show cross section vs. for the same \x < parameter space. In 

c) and d), we show the same plots, except for /i > 0. 

show points with 0.094 < h 2 < 0.129 by green dots, while blue dots denote points with 
Slzh 2 < 0.094, which would require in addition some form of non-MSSM cold dark matter. 
The corresponding plot of f-asi vs. Q^Jt? is shown in frame b.). The parameter space with 
\i > is shown in frames c.) and d.). In the case of /x < 0, much lower cross-sections are 
allowed, due to the above mentioned cancellations in quark-neutralino effective couplings. 
However, it is the /i < sign which is more restricted by and BF(b — ► S7) constraints. 

Finally we present our results on / • asi for mSUGRA together with the above men- 
tioned direct and indirect constraints in Fig. |I| for the rriQ vs. min plane. As an exam- 
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mSUGRA: tan(3=1 0, A Q =0, u>0 
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— m h LEP2 limit — a^SUSYxl 10 — Br(b^sy)x1 4 

Br(B s ^|i + n")x10 9 — 0.094<nh 2 <0.129 Stage 3 
fo^pjxio" pb — 100 .... 10 1 



mSUGRA: tan(3=45, A =0, n<0 
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m (TeV) 

m h LEP2 limit a^SUSYxl 10 Br(b->sy)x1 4 

Br(B s ^n»x1 8 0.094<ilh 2 <0. 1 29 stage 3 

fo^pjxio 11 pb 100 ....10 1 

Figure 4: Constraints and asi rates for mSUGRA 
0(upper left), tan/3 = 30,^ > 0(upper right), tan/3 = 
0(bottom right),^4 = 0. 



mSUGRA: tanp=30, A„=0, u>0 
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mSUGRA: tan(3=55, A =0, n>0 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

m (TeV) 

m h LEP2 limit a^SUSYxlO 10 Br(bn>sy)x10 4 

Br(B s ^nV)x1 8 0.094<I2h 2 <0.1 29 stage 3 

fo^pjxio^pb 100 .... 10 1 

model in mo vs. toi/2 plane: tan/3 = 10, [i > 
= 45, /.i < 0(bottom left) and tan/3 = 55, ^ > 



pie we have chosen tan/3 = 10, /z > (upper left), tan/3 = 30, \i > (upper right), 
tan/3 = 45, /x < (bottom left) and tan/3 = 55, \i > (bottom right) cases for Aq = 0. 
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The red and black shaded regions are excluded due to lack of REWSB or a stau LSP re- 
spectively. The magenta region is excluded by LEP2 searches for charginos and sleptons. 
The region below the red curve is excluded by LEP2 Higgs searches. One can see how 
BF(b — ► S7) (magenta contours), a M (blue contours), B s — > fi + fi~ (light blue contours) 
can further significantly reduce MSSM parameter space and the range of / • a si- 

For tan (3 = 10, there are three regions of allowed relic density (which could be also 
recognized by looking at different branches of Fig. |2|): 

1) lower left "bulk" region where neutralinos annihilate dominantly through t-channel 
slepton exchange; 



2) the region of t\ — Z\ co-annihilation which is along the stau LSP region [59]. One 
should notice that this region has highly fine-tuned relic density since a small variation 
in mo would lead to a large change in 42]; 



3) HB/FP region where fi becomes small and a large higgsino component of Z\ allows 



for efficient annihilation into WW, ZZ, Zh and hh pairs [61 



In the first bulk region, a si is the highest and reaches 1CP 6 — 1CP 7 pb level (Fig.||). 
It can be easily accessible by future experiments on direct CDM searches as indicated by 
yellow region in Fig. |j. Unfortunately it is already excluded by the LEP2 bound on m^. 
It also has a rather low value of BF(b — ► S7). 

The narrow fi — Z\ co-annihilation region lacks theoretical attraction due to large fine- 
tuning in the relic density, and it is also not attractive for direct DM search experiments. 
The typical level of a si is 10 -10 — 10 -11 pb (dashed and dotted black contours for tan (3 = 30 
in Fig. |4|), and in this region it will be not accessible even by Stage 3 detectors. 

Finally, from both a theoretical and experimental point of view, the HB/FP region 
looks very attractive. One can see from Fig. Q (dashed line) and Fig. |] (yellow shaded 
region) that Stage 3 experiments on direct CDM search can indeed cover the portion of 
the HB/FP region which is within the WMAP bounds. 

If we increase tan/3 to 30, as shown in Fig. [| upper-right, then most of the bulk 
region of relic density is still excluded by a combination of Higgs mass and BF(b — ► 57) 
constraints. The stau co-annihilation region remains narrow and largely beyond reach 
of direct DM detection experiments, and the HB/FP region remains consistent with all 
constraints, along with having a large rate of direct DM detection. However, by increasing 
tan j3 to 45 (// < 0) (Fig. £§ bottom-left) and tan /3 = 55 (/i > 0) (Fig. £| bottom-right) 
a qualitatively new region (fourth kind of region) of allowed relic density appears. One 
can see the diagonal strip running from the lower- left to the upper-right corner. In this 
region, neutralinos annihilate efficiently through the s-channel A and H bosons which have 
very large widths due to large b- and r- Yukawa couplings. Though these regions are also 
theoretically attractive, one can see that they will be inaccessible by neutralino direct 
search experiments. 

Finally, we would like to stress the high degree of complementarity between the LHC 
SUSY search and the direct DM search in restricting supersymmetric parameter space [20, 
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62]. A very intriguing feature of the direct DM search experiments is that they can cover 
most of the HB/FP region while LHC can explore m^ 2 only up to about 700 GeV (largely 
independent of tan/?) |63|]. On the other hand, LHC is able to completely cover the 
region of neutralino annihilation through heavy Higgs resonances and almost all of the 
stau co- annihilation region [63]. Thus, a combination of the LHC and Stage 3 direct DM 



experiments together can cover almost the entire mSUGRA parameter space! 
4. Results for Yukawa unified models 

In this Section we present results for MSSM models with non-universal boundary con- 
ditions motivated by 50(10) SUSY GUT scenarios. SUSY GUTs based on 50(10) are 
theoretically very attractive. Along with gauge coupling uification, they unify all matter 
of a single generation into a single 16 dimensional spinorial multiplet of 5O(10). The 
16 dimensional spinorial multiplet contains a right-handed neutrino which becomes a SM 
gauge singlet. It can be used to provide for massive neutrinos via the see-saw mechanism. 
The massive neutrino states may play an important role in baryogenesis via intermediate 
scale leptogenesis |34| due to the structure of the neutrino sector. Finally, 5O(10) is an 
anomaly-free group and explains the cancellation of triangle anomalies within the SM. 

In the simplest models of 5O(10) SUSY GUTs, the two MSSM Higgs doublets are 
both present in a single 10-dimensional representation. This structure of the Higgs sector 
leads to a crucial feature of the model — it predicts Yukawa coupling unification since 



the superpotential contains the term: / B fip(l6) ^(16)^(10) + • • • . In Ref. [3C], we 
have searched model parameter space for solutions with a high degree of third generation 
Yukawa coupling unification, for both positive and negative values of the superpotential \x 
parameter. In the case of \i < 0, there exist regions of SUSY parameter space satisfying 
all experimental constraints and possessing perfect Yukawa coupling unification. Yukawa 
unified solutions for fi > were also found. These latter solutions required multi-TeV values 
of the GUT scale scalar mass parameter miQ , and were in general difficult to reconcile with 
bounds on the relic density of neutralinos. 

We have adopted results from |3(]] using Isajet v7.64 and studied the neutralino-proton 
scattering cross section for cases with a high degree of Yukawa coupling unification. The 
measure of Yukawa unification used is given by 

R = max(/ (l f b , f T )/ min(/ i , f b , f T ), (4.1) 

where ft, fb and f T are the t, b and r Yukawa couplings, and R is measured at Q = Mqjjt- 
In 5O(10) models, the reduction in rank of the gauge symmetry group can lead to additional 
D-term contributions to scalar particle masses. In the case of 5O(10) — > SU(5) x U(l)x — ► 
SU (3) c x SU (2)lxU (l)y , the D-term contributions lead to the following GUT scale scalar 
mass splittings (DT model) 

m Q = m E = m U = m 16 + M D> 

r n 2 D = m L = m 16 - 3M £>> 
m 2 N = m 2 16 + 5Ml, 



-12- 



DT models with R < 1.3 and \\ < 
0<m 16 <20TeV, 0<m [J2 <5TeV, -3m <A g <3m Q , 40<tan$<60 

— 
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Figure 5: Scattering spin-independent neutralino-proton rates versus neutralino mass(left) and 
tth 2 (right) for DT(top) and iJS(bottom) models 



m 2 10 T2M 2 D , 



where parameterizes the magnitude of the D-terms. Owing to our ignorance of the 
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gauge symmetry breaking mechanism, can be taken as a free parameter, with either 
positive or negative values. \Md\ is expected to be of order the weak scale. An alternative 
scenario is for mass splittings to occur only for GUT scale scalar Higgs masses []66| (HS 
model) i.e. 



2 2 2 2 

mQ = % = mjj = m 16 , 

2 2 2 

m D = m L = m 16 , 

2 _ 2 

m N — m 16 , 
m 2 Hud = m 2 w T2Ml. 

The DT model is slightly preferred for models with [i < 0, while the HS model is preferred 
for \x > 0. Both the DT and HS models are characterized by the following free parameters, 

mi 6 , mio, Mo, m 1/2 , A , tan/?, sign(^). 

We scan these models over the following parameter range: 

< mie < 20 TeV, 
< mio < 30 TeV, 
< m 1/2 < 5 TeV, 
-(m w /V2) 2 <M 2 D < +(m 10 /v / 2) 2 , (4.2) 
40 < tan (5 < 60, 
— 3mi6 < Ao < 3mi6. 

Our results for Yukawa unified models are presented in Fig. [5] for /i < DT model 
(upper frames) and for ji > -fTS 1 model (lower frames). In each, we require R < 1.3. The 
left-most frames show / • asi vs - m z 1 ^ w hile the right-most frames show / • asi vs - ^z^ 2 - 
Green points fall within the WMAP limits on the relic density, while blue points have 
^Zi^ 2 ^ 0-094. Black points have ^ Zl h 2 > 0.129. We see immediately from the relatively 
infrequent number of green points how restrictive the WMAP bounds are on Yukawa 
unified models. Nevertheless, some Yukawa unified models are consistent with WMAP. 
These models tend to have relatively heavy sparticle mass spectra, and usually fulfill the 
£ly h 2 bound by having 2m y ~ vtla so that resonance neutralino anihilation via heavy 



and pseudoscalar Higgs bosons occurs at a high rate|3C]. The WMAP allowed models tend 
to have direct detection rates below the reach of even the Stage 3 DM search experiments. 
We do note, however, that several points do occur with allowed relic densities that are 
accessible to direct DM searches. 

In the lower two frames, we show the corresponding results for HS models with /i > 
0. In this case, values of 8 — 20 TeV are favored. The SUSY particle spectra is 
characterized by an inverted scalar mass hierarchy, so that the models may still satisfy 
naturalness criteria, since third generation scalar can be relatively light. However, the 
very heavy spectra leads usually to too high of values of fi^/i 2 , and also suppresses the 
direct detection rate to below the reach of even Stage 3 detectors. Again we note that a 
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SO(10), DT SO(10), HS 




Figure 6: Constraints and asi rates in m vs. mi/2 plane for DT(left) and iJS(right) models 

few points do survive the relic density constraint, and one is even within reach of Stage 3 
experiments. 

It is also useful to examine the Yukawa unified models in the context of constraints 
from BF(b — > s-y), and BF(B S — > fi + /i~). We illustrate this in Fig. ||, where we present 
contours of / • asi together with experimental constraints in the miQ vs. TO1/2 plane. The 
left frame of Fig. || presents results for \i < and mio = 1.356mi6, Md = 0.394mi6, Ao = 
and tan /? = 56, which leads to a central region of the plot with good Yukawa unification 
(illustrated by the R < 1.02 — 1.05 contours. One can see that the region of good Yukawa 
coupling unification can satisfy all experimental constraints with a small b — > 57 pull. The 
valid region is marginal for the future direct CDM search experiments which is indicated by 
the yellow region. The value of / • a si can be quite hight ~ 10~ 7 pb in the lower left corner 
but BF(b — ► S7) is too high in this region (magenta contour) and therefore it is excluded. 
In addition this corner of parameter space is excluded by experimental constraint on a M 
since its deviation from the SM is below = —30 x 10 -10 (blue contour) level. 

We illustrate results for the HS model with fi > in the right most frame, where 
the parameter set was chosen to give reasonable values of the neutralino relic density and 
to have Yukawa unification at the order of 20%. The left most side has a good relic 
density value but poor Yukawa unification. The low m^2 region has an area of good relic 
density (due to neutralino annihilation through the light Higgs h pole) and Yukawa coupling 
unification below R = 1.19. While much of this region is accessible to direct DM searches, 
much of it is also gives values of BF(b — * sj) outside of the window of expectation. 
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5. Gaugino mediated SUSY breaking models with non-universal gaugino 

masses 



As mentioned in the introduction, the formulation of SUSY GUTs in extra dimensions 
yields elegant solutions to many of the problems encountered by four-dimensional models [Q] . 
Recently, SUSY GUT models have been formulated in five and even higher dimensions p3[ . 
Models have been constructed wherein the doublet-triplet splitting problem is simply 
solved, and where proton decay is either suppressed or forbidden. i?-parity conservation 
may also naturally arise. All these facets may arise without the need for the large Higgs 
representations which are needed to break the GUT symmetry in 4-d models. 

In 5-d models, compactification of the extra dimension on an S 1 /(Z2 X Z 2 ) orbifold leads 
to two inequivalent fixed points identified as separated 4-dimensional branes within the 5-d 
bulk. Only fields with positive parities under the Z2 and Z' 2 symmetries have massless 
Kaluza-Klein modes. The Zi parity assignments are chosen so that the N = 1 SUSY in 
5-d, which normally reduces to N = 2 SUSY in 4-d, in fact breaks to N = 1 SUSY in 4-d. 
Superfields containing MSSM matter are assumed to exist on the observable brane, while 
the "hidden" brane is set up to accomodate SUSY breaking. The Z' 2 parity assignments 
are made so that the grand unified symmetry is broken on the hidden brane. The set up is 
well suited to accomodate SUSY breaking via gaugino mediation, wherein SUSY breaking 
is communicated from the hidden brane to the visible brane via gauge superfields which 



propagate in the bulk[67]. In gaugino mediated SUSY breaking (inoMSB), the scalar, 
trilinear and bilinear soft SUSY breaking masses are loop suppressed, and can effectively 
be taken to be zero. The gaugino masses however are non-zero. 

Several intriguing scenarios arise in gaugino mediation. In all cases, in the 4-d effective 
theory below min(M Cl Mqut), mo ~ A) ~ 0, while non-zero fj, and B can be generated. 1 
These latter parameters are traded for tan/3 and Mz, as is usual where electroweak sym- 
metry is broken radiatively. The scenarios are given as follows. 

• 1. If SU (5) is the gauge symmetry on the hidden brane, then universal gaugino masses 
will result. In this case, gaugino mediation leads to a stau LSP unless additional 



above-the-GUT-scale RGE running is allowed (Schmaltz-Skiba case) [69, 70]. 



2. If SU(5) is the gauge symmetry of the bulk, but is broken to the SM gauge 
symmetry on the hidden brane, then independent gaugino masses are induced. The 
effective 4-d theory on the visible brane is just the MSSM, with parameter space 
Mi(M c ), M 2 (M C ), M 3 (M C ), tan/3 and sign(/j,). 

3. If S'O(IO) is the bulk symmetry, and is broken to 5*0(6) x SO(<i) (isomorphic to 
the Pati-Salam group SU(4) x SU(2) L x SU(2) R ) on the hidden brane, and 50(10) 
is broken to SU (5) on the visible brane via the usual Higgs mechansim, then the 
visible sector model obeys parameter space is given by M2(M C ), Ms(M c ), tan/3 and 
sign(n), where M\ = §M 3 + §M 2 at Q = M c (Dermisek-Mafi case) [n}. 



1 Analyses of models with non-universal gaugino masses and non-zero scalar masses and ^4-terms are 
presented in Ref. |6q]. 



- 16 - 



• 4. If 50(10) is the bulk gauge symmetry, and flipped SU(5)' x U(\)' is the hidden 
brane symmetry, then the parameter space is given by Mi(M c ), M2(M C ), tan/3 and 
sign(fi), where M^{M C ) = M2(M C ) (Barr-Dorsner case) [fF^| . 

In this section, we explore the direct dark matter detection rates for the latter three cases 
of non-universal gaugino mediation. In these models, the additional parameter freedom 
arising from non-universal gaugino mediation can be exploited to solve the slepton LSP 
problem, instead of above-the-GUT scale running suggested by Schmaltz and Skiba. Hence, 
in our calculations, we adopt the choice of M c = 1 x 10 16 GeV, i.e. near but just below 
the usual unification scale of Mqut — 2 x 10 16 GeV. This choice preserves the prediction 
of gauge coupling unification. 

Our main results are shown in Fig. [j] for case 2 with independent gaugino masses at 
the compactification scale. We have scanned over the parameter space limits as listed at 
the top of the figure. Solutions which yield a spectrum with a neutralino LSP are listed by 
dots. These solutions usually have a lightest neutralino with a substantial wino or higgsino 
component, since M\{M C ) must be chosen to be large enough to contribute to RG running 
of mf so that it avoids becoming the LSP. Neutralinos with a large wino or higgsino 
component can annihilate at high rates to WW, ZZ and Zh pairs so that they generally 
give a rather low value of the relic density 1 73]. Thus, supersymmetric dark matter would 



have to be augmented by other forms of dark matter to match the constraints on ^cdm^- 2 
from WMAP and other measurements. Solutions with ^j^/i 2 > 0.129 are given by black 
dots, and would be excluded. Green dots denote solutions with 0.094 < Q^h 2 < 0.129, 
and would not require other sources of dark matter. Finally, the light and dark blue dots 
have Qg^h 2 < 0.094. The light blue solutions correspond to those with high tan/3 > 40, 
while light blue solutions have 5 < tan/3 < 40. We see from the plot that a significant 
number of solutions are accessible to Stage 1, Stage 2 and Stage 3 detectors. 

While the scans over the complete parameter space are comprehensive, it can be in- 
structive to plot results in particular parameter space planes. We present a scan of the 
M2 vs. M3 parameter space in Fig. [8| for fi > and tan/3 = 10 with M\ = 1 TeV (left 
frame) and tan/3 = 40, M\ = 1.5 TeV (right frame). It is easy to see that taking M2 
or M3 too large will result in a model with a non-neutralino LSP, while taking M3 too 
small results in a breakdown in REWSB. Taking M2 too small results in a chargino mass 
in violation of LEP2 limits. The green regions give viable solutions which have in addition 

h 2 < 0.129. The yellow regions in addition should be accessible to Stage 3 direct DM 
detection experiments. We also show various contours of values, b — > S7 branching 
fraction and a„. Whether or not any of these additional constraints are invoked, it is clear 
that at least some of the parameter space would be accessible to direct detection searches. 



If case 3 with 50(10) breaking to the Pati-Salam group is invoked, then the parameter 
space becomes more restrictive, since now M\ is determined by the choice of M2 and M3. 
In this case, we again plot parameter space in Fig. ^ in the M2 vs. M3 plane for tan/3 = 5 
(left) and 10 (right), with [i > 0. The viable parameter space is again denoted by the green 
color, with experimental constraints indicated by various contours. Much higher values 
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Figure 7: Constraints and asi rates for independent gaugino mass scenario. 



of tan/3 result in almost no viable parameter space 1 71]. The yellow regions again denote 
parameter space accessible to Stage 3 experiments. 

The M<i vs. Mi parameter space with Mi = M3 for the case of SO(10) breaking to 
flipped SU(5) is presented in Fig. 10, for tan/3 = 10 (left) and 40 (right), and for fx > 0. 
One can see that the parameter space opens up with increasing Mi, and moreover, that 
the relic density is safely below the WMAP upper limit. The region with very low values 
of M2 is excluded because the chargino mass falls below limits from LEP2. For higher 
(but still low) values of M2, the lightest neutralino is dominantly wino-like, and has a large 



- 18 - 



Independent M 1 ,M 2 ,M 3 , n>0,tan|3=10, M 1= 1 TeV 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 

M 2 (TeV) 

m h LEP2limit a p SUSYx10 10 20 60 ■ LEP2 

£2h 2 <0.129 Br(b^>SY)x10 4 2 5 Stage 3 

fo(Z,p)x10 8 pb 1 0.1 0.01 



Independent M 1 ,M 2 ,M 3 , n>0,tanp=40, M 1= 1.5 TeV 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 

M 2 (TeV) 

m h LEP2 limit a p SUSYx10 10 20 60 ■ LEP2 

■ i2h 2 <0.129 Br(b^SY)x10 4 2 5 Stage 3 

fa(Z,p)x10 8 pb 1 0.1 0.01 



Figure 8: Constraints and / • 175/ rates for independent gaugino mass scenario in [Mi — M3) 
plane 

scattering cross section with protons. However, the low value of relic density causes / • a si 
to fall below the level accessible to Stage 3 direct detection searches. 

6. Conclusions 

In this paper, we have evaluated the neutralino-proton scattering cross section in 1.) the 
mSUGRA model, 2.) in SO(10) SUSY GUT models with Yukawa coupling unification, and 
3.) in extra-dimensional SUSY GUT models with gaugino- mediated SUSY breaking and 
non-universal gaugino masses. Our results for the mSUGRA model give a comprehensive 
scan of / • a si over model parameter space using Isajet v7.65. We also compared the 
DM detection rate against other constraints including those on the neutralino relic density 
from recent WMAP data, the b — > 57 branching fraction and the muon anomalous magnetic 
moment cu. The relic density constraint yields only several viable regions of parameter 
space: the bulk region at low mo and mi /2, the stau co-annihilation region at low mo, the 
HB/FP region at large mo, and the Higgs resonance annihilation region at large tan (3. The 
bulk region contains relatively light slepton masses of 300-600 GeV, so that neutralinos 
can annihilate via t-channel graphs in the early universe. The relatively light sparticle 
mass spectra associated with this region yields observable rates for direct detection of 
dark matter, but can also lead to values of m^, BF(b — * S7) and/or in violation with 
experimental limits. The stau co-annihilation region and Higgs resonance region typically 
have heavy sparticle mass spectra, and low direct detection rates, while the HB/FP region 
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Figure 9: Constraints and / • a si rates for Pati-Salam scenario in [M2 — M3) plane 



has low \fi | values and neutralinos with a significant higgsino component. This latter 
quality allows for large neutralino annihilation cross sections in the early universe, and also 
for observable rates for direct DM detection. Direct DM search experiments can probe 
most of the viable HB/FP region. This is fortuitous, since this region is difficult to access 
via the CERN LHC or a linear e + e _ collider. Meanwhile, the CERN LHC can probe 
essentially all the Higgs annihilation region at large tan /3, and almost all of the stau co- 
annihilation region. Our analysis illustrates the complementarity of collider searches for 
supersymmetric matter with direct searches for relic DM particles. 

We also examined two classes of SUSY GUT models. The first, Yukawa unified SO(10) 
models, have non-universal scalar masses. The constraint of Yukawa coupling unification 
along with indirect constraints from relic density, b — > S7 and point to rather heavy 
SUSY spectra in these models. For \x < 0, the relic density constraint can be accomodated 
by living in the Higgs annihilation corridor. But then sparticles must be rather heavy, 
and direct detection rates are usually, but not always, low. For fi > 0, multi-TeV scalar 
masses are required, and the model is characterized by a radiatively induced inverted mass 
hierarchy. The very heavy scalars suppress neutralino annihilation cross sections, so the 
relic density is typically very large, and direct detection rates are very low. 

The other class of models, extra dimensional SUSY GUTs with non-universal gaugino 
masses, require mo ~ Aq ~ 0. But then the compactification scale gaugino mass M\ 
must be typically larger than M2 and M3 to avoid a charged LSP. In this case, the LSP 
can be a neutralino, but with a significant wino or higgsino component (in mSUGRA, the 
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Figure 10: Constraints and / • asi rates for nipped SU(5) model in (M 2 — M 3 ) plane 



LSP is almost always bino-like). This leads to large neutralino annihilation rates and a 
low relic density, but it also leads to large direct detection cross sections. A portion of 
the parameter space of the independent gaugino mass case and the Pati-Salam case can 
be probed by direct detection experiments. However, in the case of flipped SU(5) gauge 
symmetry on the hidden brane, the relic density is so low that the rescaled direct detection 
cross sections are almost always below the level of even Stage 3 experiments. 
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